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ABSTRACT 

We present results from a total of 459 repeated 3.1 GHz radio continuum observations (of which 379 
were used in a search for transient sources) of the ELAIS-N1, Coma, Lockman Hole, and NOAO Deep 
Wide Field Survey fields as part of the Pi GHz Sky Survey (PiGSS). The observations were taken 
approximately once per day between 2009 May and 2011 April. Each image covers 11.8 square degrees 
and has 100" FWHM resolution. Deep images for each of the four fields have rms noise between 180 
and 310 /iJy, and the corresponding catalogs contain ~ 200 sources in each field. Typically 40 - 50 
of these sources are detected in each single-epoch image. This represents one of the shortest cadence, 
largest area, multi-epoch surveys undertaken at these frequencies. 

We compare the catalogs generated from the combined images to those from individual epochs, 
and from monthly averages, as well as to legacy surveys. We undertake a search for transients, with 
particular emphasis on excluding false positive sources. We find no confirmed transients, defined here 
as sources that can be shown to have varied by at least a factor 10. However, we find one source which 
brightened in a single-epoch image to at least six times the upper limit from the corresponding deep 
image. We also find a source associated with a z = 0.6 quasar which appears to have brightened by a 
factor ~ 3 in one of our deep images, when compared to catalogs from legacy surveys. 

We place new upper limits on the number of transients brighter than lOmJy: fewer than 0.08 tran- 
sients deg -2 with characteristic timescales of months to years; fewer than 0.02 deg -2 with timescales 
of months; and fewer than 0.009 deg -2 with timescales of days. We also plot upper limits as a function 
of flux density for transients on the same timescales. 

Subject headings: catalogs — methods: data analysis — radio continuum: galaxies — radio continuum: 
general — surveys 



1. INTRODUCTION 
1.1. Recent Results from Radio Transient Surveys 

The past few years have seen an increase in multi-epoch 
surveys specifically designed to probe the population of 
radio transients, as well as transient searches under taken 
using archiva l data (for a recent summary, see | Fender | 
& Bell 2011[ ). Transients can arise f rom processes tak - 
ing place m a variety of progenitors ( |Lazio et al.| [2009), 
with typical timescales dictated by the characteristic pro- 
cesses causing the transient emission, and the sizes of the 
emitting regions. Radio transients are generally classi- 
fied as short (characteristic timescales < 1 s) and long 
(timescales > 1 s) , where the shorter timescales typically 
correspond to coherent emission pro cesses and the longe r 
timescales to incoherent emission ( |Cordes et al.| p004). 
The timescales probed correspond to characteristic pro- 
cesses and size scales in the emitting objects, and so two 
surveys with differing temporal coverage will see a differ- 
ent mix of objects from each other. 

Of particular interest for the current generation of ra- 
dio sky surveys are tidal disruption events (TDEs) — 
flares connected to sudden increases in accretion rate 
caused when stars are tidally shredded as they pa ss close 
to black holes with mass < 10 8 M© ( |Rees||l988| . Such 
events may be the dominant contributor to transient 
source counts in th e radio for surveys focusing on long 
duration transients ( |Frail et al.||2012| Bower|[2011[ ). Re- 



cently axomioM^OTDE candidates have been detected 
(e. g. , |Gezari et aL)2009| |van Velzen et al.||2011| ) by high 
energy surveys, ol w hich some have been" detected in 
radio -follow-up (e. g. , Zauderer et al.pOlT Cenko et al. 



2012 ). However, no TDEs have yet been unambiguously 
detected in blind radio surveys. 

Other classes of sources that are expected to dominate 
blind transient surveys are radio supernovae, gamma 
ray burst afterglows, radio counterparts to gravitational 
wave sources such as neutron star mergers, flares from 
low-mass stars and brown dwarfs, and refraction and 
scintillation of extragalactic sources due to dumpiness 
and turbulence in our own Galaxy's interstellar medium 
(|Frail et al.|2012| |Bower et al.||2011|) . Some of the afore- 



mentioned source classes have been detected in follow-up 
of triggers at other wavelengths, but until recently, radio 
surveys have not had the required throughput to detect 
them blindly in significant numbers. Some progenitors 
may only be detectable in the radio, however, and in 
other cases radio variability may be more extreme than 
variability at other wavelengths, making these sources 
easier to find. Blind surveys also have the potential to 
serendipitously detect new classes of objects. 

Some transient candidates have been detected in blind 
radio surveys , thou gh — for example, those reported by 
Bower et al. ( |2007[ ), a larg e fraction of which remained 



undetected in dee p optical (|Bower et al.|2007 Ofek et al. 
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2010| and X-ray ( Croft et al.|2Ullb ]) follow-up. A recen t 
reanalysis of the Bower et al.]data py Frail et al. (2012) 



however, suggests that some of these transients may be 
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data artifacts, and the significance of s ome of the others 
may have been overestimated (Section [Q| ). Many other 
radio transient surveys may also have been somewhat 
overoptimistic in their claims of detections, particularly 
of sources seen in a single epoch and never again, and 
we urge appropriate skepticism. In this paper we discuss 
some of the safeguards we have implemented to ensure 
that any transients we detect have a high probability of 
being real. 

1.2. Limitations of Current Surveys 

Survey design involves trade-offs between a number of 
factors, including frequency, bandwidth, frequency res- 
olution, spatial resolution, limiting sensitivity, sky area, 
number of epochs, cadence (the time between epochs), 
integration time at each epoch, and Galactic latitude. 
The number of transients detected by a survey depends 
on the lightcurves, radio spectra, luminosity functions, 
and distribution in Galactic latitude and in distance, for 
a mix of different classes of progenitor, so comparing re- 
sults from one survey to another is not trivial. 

For example, a survey with many epochs with a ca- 
dence of days might see different objects to a shallow 
large-area survey with fewer epochs with a timescale of 
months, even if the effective area covered (the number of 
pointings, multiplied by the area covered by each point- 
ing) is similar. Large area surveys with contiguous cover- 
age have the advantage that sensitivity is approximately 
uniform across most of the resulting mosaic, whereas sur- 
veys with small numbers of pointings are more affected 
by primary beam falloff at the mosaic edges. Small area 
surveys with many repeated epochs, however, have the 
advantage of allowing a deep comparison image to be 
constructed from all of the data. Different authors also 
define "transient" in different ways, and as a result, mea- 
surements of transient surface densities (sometimes re- 
ferred to as "two-epoch rates" of transients per unit ef- 
fective area) from different surveys may not be directly 
comparable. 

The sensitivity of a given survey is particularly criti- 
cal in influencing the number of transients seen; a tran- 
sient seen in one survey may be classified as a variable 
source by a more sensitive survey if detected in additional 
epochs. 

We suggest that a sensible threshold for a source to be 
considered a transient is that it brightens by a factor of at 
least ten. The distinction between transient and variable 
sources is always somewhat arbitrary, but we pick this 
threshold to distinguish between sources that would be 
detected in quiescence with a modest increase in sensitiv- 
ity (e.g., variable AGNs) and those that brighten more 
dramatically due to some explosive event (e.g. SNe), 
very significant changes in accretion rate (e. g. , TDEs), 
and the like. As shown in |Bower et al.| (2011), and in 



Section |5.1[ a few percent ol sources exhibit apparent 
variability of a factor of a few, but almost none of these 
"normal" variables exhibit variations of an order of mag- 
nitude or more in flux. The choice of a factor 10 for 
transient detection naturally limits the sensitivity of any 
two-epoch survey to around 30<r (for a source to be un- 
detected at < 3<j in one image, and to change in flux 
density by a factor of 10 or more on its detection in the 
other image). 

It should also be noted that explosive events may oc- 



cur within galaxies that are themselves radio sources (for 
example, radio galaxies). If the spatial resolution is not 
sufficient to distinguish between the two sources of ra- 
dio emission, then the explosive event will not be char- 
acterized as a transient, by our definition, unless it is 
at least nine times brighter than the emission from the 
host galaxy. This will also tend to reduce the amount of 
variability seen by surveys like ours (with relatively low 
spatial resolution) where extended non- variable emission 
is unresolve d from variable so urces such as radio galaxy 
cores (e. g. , |Duric et al.|1989| ). 

The distinction between transient and variable sources 
aside (which is of most importance in comparing tran- 
sient statistics between surveys) it is of course still of 
interest to study sources which exhibit large changes in 
flux density, and we remark upon some of these extreme 
variables, as well as analyzing variability statistics for 
large numbers of sources, in Section [5] 

If we consider a multi-epoch survey with a detec- 
tion threshold of 6cr S i n gi e for a transient in an individ- 
ual epoch, and a Zodeep upper limit in a deep com- 
parison image made from all observations of the same 
field (with a deep — °~ single / \fn f° r n epochs), we see 
that to detect transients close to the threshold (where 
6a 'single > 10 x 3a 'deep) requires 25 or more epochs. Such 
repeated observations of the same field (as opposed to 
using the allocated observing time to increase the area 
of the survey while maintaining shallower depth) also 
have the advantage of allowing us to extract multi-epoch 
lightcurves for sources of interest, and to make images 
from subsets of the data with different timescales. Re- 
peated observations of the same field also allow us to dis- 
entangle intrinsic variability from telescope calibration 
errors. However, these requirements must be balanced 
against the fact that quadrupling the area covered by a 
survey at the expense of doubling the detection threshold 
should result in an increase of the number of transients 
detected of ~ 40% (naively assuming an S~ 3 ^ 2 depen- 
dence of rates on flux density, S) . The tradeoff between 
depth and area must also be balanced with the desire to 
obtain adequate sam pling for the timescales of interest. 

Some surveys (e.g. |Croft et al.||2010 ) use a deep com- 
parison field that was taken a long time (compared to the 
timescales of interest) before the survey data. For sur- 
veys which use their own data to generate a deep field, 
there is an additional complication. Transients that are 
bright, or those that are detected in more than a few 
of the individual images, will also contribute significant 
flux to the deep image, potentially resulting in a detec- 
tion even when no long-term counterpart is present. 

1.3. False Positives 
The reassessment o f the iBower et al. l (120071) transients 



Frail et 



all (120121); and of th e Levinso n^aT] (|2002[) 
By iGal-Yam et al.l (|2UUbb a nd |Ulek et al. 
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transie nts 

J2010 ); as w ell as questions raised by Burke- Spolaor 
let al.| ([20 11 1 about transients such as those seen by' 



Lori mereTaL] ( |2007 ) 
( |2009p transientTBy 



and about the |Matsumura et ah] 
Croft et al.| ( |20lia| ), illustrate that 
events detected only in a single epoch have a significant 
chance of being spurious. In fact many current transient 
searches may be operating in an uncomfortable regime 
where the number of false positive detections is larger 
than the number of true transients. 
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Spurious sources in imaging searches can arise due to 
amplification of noise by overly large primary beam gain 
corrections; from various image artifacts including side- 
lobes of bright sources and corrupted data; or simply 
due to thresholds that are overly optimistic. The latter 
can occur because the large area surveyed means that 
even rare events are seen occasionally; because underes- 
timates of the background noise cause the significance of 
detections to be overestimated; or because the assump- 
tion of Gaussian noise is not valid (which it may not be 
in the presence of artifacts as well as non-linear imag- 
ing algorithms such as CLEAN). If the primary beam is 
well-known, relatively conservative gain cutoffs are easy 
to implement; in this paper we discard sources where 
the pri mary beam correction is greater than a factor 5 
(Section |3.3 ). 

To have the best chance of finding real transients, we 
must carefully choose thresholds to lower the number of 
false positives. One way is to choose a very conserva- 
tive threshold, perhaps 8 - 10a, in individual epoch im- 
ages (and look for counterparts at any flux density in the 
deep image). However, the assumed S~ 3 ^ 2 dependence 
of rates on flux density, 5, means that high thresholds 
also sacrifice a lot of sensiti vity to potential transients. 

A better method (Section 3.5 ) is to split the (ix, v) data 
in half and create two separate images, which are then 
searched simultaneously for candidates appearing in the 
same place. Noise peaks will tend not to coincide in two 
images, and since sidelobes and other image defects also 
often do not coincide, this method can also reduce the 
number of false positives due to artifacts. 

The data can be halved in a variety of ways — in time, 
frequency, or by taking subsets of baselines, for exam- 
ple. Then the catalog can be generated with a lower 
threshold in both images (~ 4.2<j, in our case). This 
will result in relatively large numbers of spurious sources, 
but any source not appearing in both images is then dis- 
carded, resulting in a comparatively clean catalog with a 
~ 6<j threshold. For the number of independent beams 
searched, the probability of a single source appearing in 
our catalogs above this threshold due to random Gaus- 
sian noise is 0.7%; as noted above, however, artifacts or 
other imaging defects can still easily result in the pres- 
ence of spurious detections above this threshold. 

The observing strategy used for PiGSS (two correlators 
delivering adjacent frequency bands) lends itself natu- 
rally to the two-image approach, but it could be applied 
by other telescopes with a single correlator too, by di- 
viding the data up in some way before imaging. Divid- 
ing in frequency may of course miss transients with very 
steep spectral indices, and dividing in time reduces sen- 
sitivity to rapid variability, but in such cases more than 
one method of halving the data could be used, and the 
catalogs checked against each other. We discuss our im- 



plementation of this method in Section 3.5 



2. TRANSIENT SURVEYS WITH THE ALLEN TELESCOPE 

ARRAY 



The Allen Telescope Array (ATA; [Welch et al.| [2009), 
is a 42-element interferometer, located at Hat Creek Ka- 
dio Observatory in Northern California. Each element 
is a 6.1-m dish with a broadband feed that receives ra- 
diation from 0.5 to 11 GHz. Two independently-tunable 
correlators select bands of 104 MHz width, each consist- 



ing of 1024 channels, from anywhere in this frequency 
range. The maximum baseline is 300 m, yielding resolu- 
tion ~ 100" at ~ 3 GHz. Typical system equivalent flux 
densities at ~ 3 GHz for antennas that are performing as 
e xpected are y 10 4 Jy. 

Croft et al.M2010|) andlCroft et al.M2011a|) presented re- 



sults from a pilot program, the ATA Twenty-centimeter 
Survey (ATATS). They compared the resulting cat alogs 
to those from th e NRAO VLA Sky Survey (NVSS; [Con^ 



don et al. 
with the 



1998), a survey also at 1.4 GHz undertaken 
VETTbetween 1993 and 1997. IBower et a~ 



(2010), hereafter Paper I, introduced a larger ATA sur- 
vey, the Pi GHz Sky Survey (PiGSS). PiGSS covered 
~ 2500 square degrees in two epochs, and an additional 
~ 2500 square degrees in a single epoch, in two simul- 
taneous 100 MHz bands centered at 3.04 and 3.14 GHz 
(Croft et al. in prep.). As well as the large field, PiGSS 
also incorporates four smaller (~ 12 square degree) fields 
observed with a cadence of days. Paper I presents results 
from a deep image made from 75 epochs of observation of 
one of these fields (centered on the l ocation of the NOAO 
Deep Wide Field Survey, NDWFS). [Bower et al | ( [20TT] ), 
hereafter Paper II, compares catalogs generated from the 
individual epochs of this field (and mont hly averages of 



them ) to each other, and to other radio (iCondon et al. 
" jBecker et al.|[l995| [Gregory et al.||1996| |de Vries 

Kc 



infrared (Ashby et al. 2009[" Kozlowski et al. 
K-ray (Kenter et al. 2005 p catalogs. 



et al. | 2002 | ) 
[201% and X 

In this paper, we present results from the other three 
12 square degree deep fields, centered on the ELAIS-N1, 
Lockman Hole, and Coma fields (Table [I]), as well as a 
partial reanalysis of the NDWFS results from Paper II 
to conform with updated analysis techniques used on the 
other three fields. The ELAIS-N1 and Lockman fields 
were selected due to the availability of large amounts 
of deep imaging and spectroscopic coverage from other 
telescopes. The Coma field also has extensive multi- 
wavelength coverage, and in addition the wide- field of 
the ATA enabled our survey to include a large number 
of low-redshift galaxies in the Coma cluster, as opposed 
to the higher-redshift galaxies in the other fields. All four 
fields are at high Galactic latitude (b > 50°) and as such 
primarily target extragalactic transients; results from an- 
other recent ATA survey primarily targeting G alactic 
transients are reported by Williams et al. (2012). Most 
existing transient surveys have to cussed on frequencies in 
the 1-2 GHz range, with typical cadences rather longer 
than used for our observations, and so PiGSS provides a 
unique view of the transient and variable radio sky. 

We place a particular emphasis here on techniques that 
will be scalable to the next generation of imaging surveys, 
where processing and analysis must take place with an 
increasing degree of automation due to the very large 
data volumes involved. Although PiGSS is still (just) at 
a scale where manual processing is possible, we focused 
on developing routines to reject bad data both before 
and after the imaging stage, and on understanding our 
results in the context of other surveys, both past and 
future. 

Throughout this paper, we use J2000 coordinates. 
3. DATA ACQUISITION AND REDUCTION 

The PiGSS observations reported here were under- 
taken approximately daily between 2009 December and 
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Table 1 

PiGSS Deep Fields 



Field 


RA 


Dec 


Date Range 




Number of Epochs 




Deep Field 




(J2000) 


(J2000) 




Total 


Deep a 


Transient 


Variable 


rms (/xJy) 


NDWFS d 


14 32 00 


34 16 00 


2009 May 20 Sep 29 


96 


78 


67 


66 


250 


Lockman 


10 52 00 


57 38 00 


2009 Dec 2 2010 Jul 7 


132 


87 


110 


101 


250 


ELAIS-N1 


16 12 42 


53 48 31 


2010 May 30 2011 Jan 20 


179 


154 


165 


145 


180 


Coma 


13 00 43 


27 11 13 


2011 Jan 24 Apr 7 


52 


43 


37 


33 


310 



a Median number of epochs for each field with rms flux density < 3.5 mJy beam -1 in images made from individual 
pointings. These data were used to generate the deep field images, except for the NDWFS deep field where selection was 
performed differently (see Paper I) 

b Number of epochs with 90% completeness better than 50 mJy at both frequencies, that were used in the transient search 
(Section [4}. 

c Number of epochs with 90% completeness bett er th an 50 mJy at both frequencies, and Post Imaging Calibration gain 
within 0.1 of the median for that field (Section |3.4| ). These were the epochs used to determine variability statistics 
(Section [5]). 

d Reanalysis of data presented in Paper I and Paper II. 



2011 April, in campaigns lasting several months for each 
field (Table [lj. We attempted to observe all seven point- 
ings for each field for three minutes, three times, dis- 
tributed in hour angle so as to optimize (ii, v) coverage, 
during the 4 - 12 hr allocated to each PiGSS epoch. As 
in Paper I, pointings were separated by 47', the volt- 
age FWHM at 3.14 GHz, to provide approximately uni- 
form sensitivity in the central region of the field. Data 
were taken in two 100 MHz bands centered at 3.04 and 
3.14 GHz, each covering 1024 channels. One of the bright 
calibrators 3C 286, 3C 48, or 3C 147 was observed hourly. 

3.1. Flagging 

Bad data arose from a variety of sources: radio fre- 
quency interference (RFI); clock synchronization prob- 
lems; Walsh system problems; temperature related prob- 
lems for individual boards (particularly radio frequency 
converter boards) causing spectral corruption; antennas 
that were offline or had problems with feeds; and other 
issues, known and unknown. Most problems except RFI 
were relatively infrequent and typically affected a single 
epoch of data. RFI occurred in essentially all datasets, 
both as narrow-band spikes (relatively easy to excise) and 
wide-band RFI affecting a larger number of frequency 
channels. Identifying and excising ("flagging") these bad 
data in an efficient, automatic manner is important for 
large surveys like ours where the large amount of data 
makes pipeline proce ssing essential. F lagging was per- 
formed by MIRIAD flSault et aLl|1995| ) tasks wrapped 
in our own custom Perl pipeline ( autoflag") which ex- 
pands and improves upon the earlier RAPID / ARTIS 
( Keating et aLpOlO ) ATA data reduction software. 

The autoflag software computes the median and vari- 
ance of the real and imaginary values for each channel / 
polarization / baseline combination in each single-epoch 
dataset as a function of time. Channels where more than 
60% of the visibilities deviate by more than four standard 
deviations from the median values, as well as their im- 
mediate neighboring channels, are flagged as suspected 
RFI. Additionally, individual visibilities that are more 
than 10cr from the median amplitude are flagged. To 
flag wide-band RFI, spectra are binned (using a range 
of bin sizes) and bins with excessive variance (three or 
more times the median) are flagged, on a baseline by 
baseline basis. Where exactly one quarter of the chan- 
nels in spectra associated with a particular antenna show 



large variance, spectral corruption is suspected and the 
channels are flagged on all baselines associated with that 
antenna. Baselines are also flagged if calibrator phases 
show excessive (> 1 radian) variability. 

Baselines shorter than 110 wavelengths (which tended 
to be subject to solar interference), along with 100 chan- 
nels at each edge of the band (due to filter roll-off) and 
the central DC channel were also flagged. The visibility 
retention (the percentage of cross-correlations for each 
dataset remaining after all flagging) ranged from 1% to 
78%. 93% of datasets had retention of at least 65%. 
Fields with poor retention often showed obvious defects 
after imaging, or failed to image entirely; such fields were 
excluded from further analysis. 

3.2. Calibration and Imaging 

Where multiple calibrators were used, the primary cal- 
ibrator was taken to be the one with the most unflagged 
data, and the reference antenna was taken to be the one 
with the most unflagged visibilities in the primary cali- 
brator data. Initial calibration was performed using MF- 
CAL with flux densities from |Baars et al.| ( |1977| ), and 
then amplitudes which deviated by more than 4 Jy from 
the nominal flux density, and phases which deviated by 
more than 20°, were flagged in the calibrator data. A 
second round of MFCAL was performed, bandpass so- 
lutions were copied to the other calibrators (if present), 
and MFCAL was run on these calibrators to obtain a gain 
solution. The gain and bandpass solutions were copied 
to the target data from the primary calibrator, and the 
gain solutions for the other calibrators (if present) were 
merged into the target data. No polarization calibration 
was performed. 

At each pointing, images of 1024 x 1024 pixels, where 
each pixel is 15 x 15", were produced, using natural 
weighting and multi- frequency synthesis. In order to pro- 
duce mosaics, the 'mosaic' option was used with the IN- 
VERT command, and the 'offset' keyword used to set 
the reference coordinate for all seven pointings to a sin- 
gle position at the center of the mosaic. If these steps 
are not performed, MIRIAD LINMOS does not correctly 
handle the sky projection when mosaics are produced, 
resulting in offsets between ATA positions and the stan- 
dard reference frame. This was the cause of the offsets 
between ATA and NVSS positions seen in Paper I, which 
are corrected here in our reanalysis. 
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An initial CLEAN was performed down to the first 
negative component, and sources in the resulting map 
were found using the MIRIAD source finder, SFIND. The 
output image from SFIND was convolved with the dirty 
beam, and pixels below 10<j were blanked. This image 
was then used as a mask for a second round of CLEAN 
(on the image which had already been CLEANed down 
to the first negative component), which continued un- 
til a cutoff of 5mJy or 5000 iterations was reached. 
The two stage CLEAN ameliorates non-linear behavior 
in CLEAN, and essentially performs more iterations on 
sources brighter than 10<j in the initial image. 

If sources brighter than 500 mJy were present in the 
field, amplitude self-calibration was performed, and then 
the image was inverted and CLEANed again. The final 
image was restored using a 100" x 100" Gaussian beam. 

This reduction process generated over 6000 images — 
two (one for each frequency band) for each of the seven 
mosaic pointings at each of the tens of epochs for each 
of the four fields. We computed the background rms 
for each of these images using a skewed Gaussian fit to 
the pixel histogram, with iterative clipping to remove 
source pixels. Those of the ~ 6000 datasets which re- 
sulted in images with rms < 3.5 mJy beam -1 (Fig. [2| 
were considered to be good datasets, for the purposes 
of creating deep images (Figs. [4j [5j and [6] except that 
the NDWFS deep image was constructed using the same 
single-epoch datasets used in Paper I). One image was 
made for each pointing for each frequency in each field. 
Due to small variations in image quality, the number of 
epochs which satisfied the 3.5 mJy beam -1 threshold var- 
ied by ~ ±10% across the seven pointings, resulting in 
slight (~ 5%) variations in limiting sensitivity across the 
mosaic, which we neglect in subsequent analysis. In Ta- 
ble [l] we report the median number of epochs that went 
into each pointing. 

The resulting seven all-epoch single-pointing images 
at each frequency were linearly combined into two all- 
epoch single frequency band "master mosaics", using a 
tapered primary beam correction in LINMOS. We also 
linearly combined these pairs of mosaic images into a 
dual-frequency mosaic for each field. 

In addition to these three deep images for each of the 
fields, we made seven-pointing mosaic images of each 
epoch at each frequency, and linearly combined these into 
a dual-frequency mosaic for each field. Typical images, 
here de fined as mosaics with the median completeness 
(Section 3.6) for each field, are shown in Figs. HI |5l and 



[6] We also made seven-pointing mosaics of all goocTdata 
taken in a given month for each field, again inverting 
data for each pointing separately, and mosaicking in the 
image domain with a tapered primary beam. 

Reprojection to an equal- area coordinate system is nec- 
essary for large mosaics, in order that flux densities are 
measured correctly. However, for the relatively small 
fields studied here this effect is negligible (< 0.1% varia- 
tion in pixel area in the default SIN projection), and we 
chose not to reproject in order to avoid possible aliasing 
as well as the additional processing required. 

The single-epoch, single-pointing images had a median 
rms of 2.42 mJy (compared to theoretical expectations of 
~ 1.4 mJy). The single-epoch mosaics had a median rms 
of 2.02 mJy. The monthly images had a median rms of 
0.43 mJy. The deep images had rms of a few hundred 



/iJy (as expected when scaling the rms of the individual 
images by the square root of the number of epochs), as 
shown in Table [T] The ELAIS-N1 field reaches a sensitiv- 
ity close to our estimate (Paper II) of the rms confusion 
limit, ~ 150 /iJy. 

Image quality can be assessed in a variety of ways. We 
inspected all of the images by hand to ensure that select- 
ing particular methods and thresholds did not inadver- 
tently result in the inclusion of bad data (or the exclusion 
of significant amounts of good data). For the next gener- 
ation of large surveys, however, inspecting more than a 
small fraction of the images will become impractical. For- 
tunately, good images have certain characteristics: they 
have noise that is reasonably low, they contain a reason- 
able number of real sources (Fig.Jl]), and the sources have 
reasonable flux densities. 

As discussed above, we generated deep field images 
based on the noise measured in images made from data 
for individual pointings, keeping only images with rms 
noise < 3.5 mJy beam -1 . As shown in Fig.|2j and in Ta- 
ble [TJ around 85% of the ELAIS and Coma single-epoch 
images satisfy this criterion, and around 65% of the Lock- 
man images. Some epochs have rms higher than this cut- 
off, but can nevertheless be used to search for transients 
and variable sources, whereas a handful of epochs are 
completely unusable. 

To determine transient source counts or limits we wish 
to determine the effective area covered by the survey as 
a function of completeness limit (Section 3.6). Images 
which appear cosmetically very poor (some ol which also 
contain spurious sources in the catalogs as a result of 
image defects) will also tend to have poor completeness 
limits. For our transient search, we discard epochs where 
90% completeness is not reached by 50mJy (Fig. [3]). The 
number of epochs surviving these cuts for each field are 
shown in the "Transients" column in Table [TJ 

To measure variability, we require relatively good com- 
pleteness for all images, as well as reliably measured flux 
densities. For our variability analysis, we adopt a com- 
pleteness threshold of 50mJy, and r equi re post imaging 
calibration gain corrections (Section 3.4) within ±0.1 of 
the median value for that field. The number of epochs 
surviving these cuts for each field is shown in the "Vari- 
ability" column in Table [I] 

3.3. Catalog Creation and Matching 

Source detection, matching, catalog creation, source 
databases and matching, and lightcurve and postage 
stamp creation were performed by our custom software 
known as "SLOW", a suite of Python and Perl scripts 
calling MIRIAD routines and building and querying 
sqlite3 databases. 

SLOW first runs MIRIAD SFIND on the master mo- 
saic and on the single-epoch mosaics. W e used SFIND' s 
false detection rate (FDR) algorithm (Hopkins et al. 
2002 ) with a box size of 50 pixels, and the default 2% 



acceptable percentage of false pixels. This corresponds 
to a detection threshold of around 4.2a. The "psfsize" 
option was used, which restricts the minimum fitted size 
of sources to the size of the restoring beam; without this 
option, some sources in the single-epoch images tend to 
be fit by ellipses with minor axis < 100", causing their 
flux densities to be underestimated relative to those from 
the master mosaic. Sources poorly fit by SFIND (indi- 
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Figure 1. Number of sources detected by SFIND in images of all individual single-epoch mosaics at 3040 (solid line) and 3140 MHz 
(dashed line) for ELAIS-N1, Lockman, Coma, and NDWFS. 
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Figure 2. Flux density rms for images of all individual pointings at 3040 (solid line) and 3140 MHz (dashed line) for ELAIS-N1, Lockman, 
Coma, and NDWFS. The dotted lines show the 3.5 mJy beam -1 cut used to select data to construct the deep images. 
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Figure 3. 90% completeness limits for all individual single-epoch mosaics at 3040 (solid line) and 3140 MHz (dashed line) for ELAIS-N1, 
Lockman, Coma, and NDWFS. The dotted lines show the 50mJy cut used to select epochs for the transient search. 
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Figure 4. Left: Lockman deep field, made with data at both 3040 and 3140 MHz. The greyscale runs from to 10 mJy beam 1 . Right 
Typical epoch (90% completeness 15.6 mJy) of the Lockman field at 3040 MHz. The greyscale runs from to 50mJy beam -1 . 
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Figure 5. Left: ELAIS-N1 deep field, made with data at both 3040 and 3140 MHz. The greyscale runs from to lOmJy beam 1 . Right 
Typical epoch (90% completeness 10.0 mJy) of the ELAIS-N1 field at 3040 MHz. The greyscale runs from to 50mJy beam -1 . 
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Figure 6. Left: Coma deep field, made with data at both 3040 and 3140 MHz. The greyscale runs from to lOmJy beam 1 . 
Typical epoch (90% completeness 15.2 mjy) of the Coma field at 3040 MHz. The greyscale runs from to 50 mJy beam -1 . 
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cated by asterisks in the SFIND output), typically com- 
prising ~ 3% of the total sources detected in each epoch, 
were rejected from the catalog. We computed total flux 
density uncertainties by adding the reported background 
rms and fit uncertainties from SFIND in quadrature. 

We corrected the attenuation at the edge of the mosaics 
by creating a gain map for each epoch in LINMOS, and 
correcting the measured flux densities and uncertainties 
in the catalogs for the appropriate gains (since SFIND's 
built-in primary beam correction does not work for linear 
mosaics). Sources in regions with gain corrections greater 
than a factor 5 were discarded from the catalog. This 
helps to minimize the problem of noise variations close to 
image edges with large gain corrections being interpreted 
as real sources. 

Sources were added to the database one epoch at a 
time, and paired with existing sources if there was a 
match within 150". The full sqlite3 databases are avail- 
able from the first author on request. 

3.4. Post-imaging calibration 

We performed a modified version of the p ost-imaging 
calibration (PIC) of Bannister et al. (2011) on our cat- 
alogs. Sources in the individual epoch catalogs were 
matched to those from the master catalog, and a linear 
least squares fit to the primary beam corrected flux den- 
sities was performed. An ordinary least squares fit (with 
weights given by the image uncertainties) gives good re- 
sults for the majority of epochs, but we were concerned 
that real variations in source flux densities could cause 
some epochs to have poor fits. While this should not be 
a problem for transient detection (since transients will 
not have a match in the master catalog), sources with 
extreme flux density variations could cause poor fits. To 
test this, we made a copy of the SLOW database which 
we then manipulated for testing purposes. For a source 
detected in most of the individual epochs, and also de- 
tected in the master mosaic with flux density 15mJy, 
we increased its cataloged flux value at each individual 



epoch to 1000 mJy, and then ran the fitting procedure 
again. Figure [7] shows data from a typical epoch, where 
this artificially modified source (designed to simulate an 
extreme variable or a source with an erroneous flux den- 
sity measurement in real data) causes the least squares 
fit to fail. Performing a RANSAC fit flFischler fc Bolles 
1981), however, results in fits to the data which are very 
robust to outliers, and so we adopt this technique for all 
of our fits. No sources in our data show variations in flux 
density as extreme as our simulated example, so our fits 
should describe the data very well. 

If we fit for both the slope and intercept, we find that 
the absolute value of the median intercepts are < 2 mJy 
for the Lockman and ELAIS-N1 fields at both frequen- 
cies, which suggests that CLEAN bias is not significant. 
For the Coma field, median intercepts are ~ — 6mJy. 
Coma exhibits extended structure and a bright source, 
which make imaging and source detection more challeng- 
ing and probably contribute to the worse fits and possible 
CLEAN bias here. For subsequent analysis, we follow 



Bannister et al. (2011), and constrain the fits to pass 



through zero. 

We use the slope of the least squares fit to estimate 
a PIC gain correction for the single epoch mosaics rela- 
tive to the master mosaic. We obtain median PIC gains 
of 0.94 for Lockman, 0.98 for ELAIS-N1, 0.96 for ND- 
WFS, and 0.89 for Coma. PIC gains smaller than one 
mean that single-epoch flux densities somewhat under- 
estimate deep image flux densities. The standard de- 
viation (measured for each field and frequency) of PIC 
gains is ~ 0.1 (e.g. , Fig. [8|, suggesting that the overall 
calibration is good to about this level. We can also use 
these gain variations to estimate a cutoff for the convinc- 
ing detection of variability in lightcurves; the epoch-to- 
epoch variations in gain suggest that a sensible threshold 
(for well-detected sources) is 3a ~ 30%, although some 
sources do have discrepant measurements that are larger 
than this, even in epochs where the PIC gain suggests 
an overall calibration that is reasonable. In a small frac- 
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Figure 7. Flux densities of sources detected in both a typical 
epoch (epoch 84) and also in the master catalog of PiGSS ELAIS- 
Nl 3040 MHz. One of the sources, detected with 19.4 mJy in the 
master catalog, has had its flux density in the single epoch catalog 
manipulated to make it into an outlier with 1000 mJy, in order to 
test the robustness of our fitting procedure. The dashed line shows 
the results of a weighted linear least squares fit: a slope of 0.75, and 
an intercept of 33 mJy, which were influenced by this single outlier. 
The solid line shows the RANSAC fit: slope 1.03, and intercept - 
2mJy, which are more representative of the inliers in the dataset. 
When determining fits for the entire dataset, we use the RANSAC 
method with fits constrained to pass through zero. 

tion (~ 7%) of epochs with a large amount of scatter 
between the single-epoch and deep-image flux densities, 
the RANSAC fit fails to converge to within our chosen 
parameters (at least 70% of sources selected as inliers af- 
ter 1000 iterations) and these epochs are rejected in our 
variability analysis (Section [5]). 

Since the PIC gain corrections are relatively small (al- 
beit larger in Coma due to the extended structure and 
bright sources), and uncertain themselves to at least a 
few percent due to the small number of sources present 
in the fits, the intrinsic variability of individual sources, 
and the variation in thresholds from epoch to epoch, we 
do not apply corrections to fluxes measured in individ- 
ual epochs. The PIC gains do, however, provide a good 
means by which to identify bad epochs. For variability 
studies, PIC gain stability is relatively important, lest 
PIC gain variations be mistaken for intrinsic source vari- 
ability (Section^. For transient searches, larger discrep- 
ancies in PIC gams may be tolerated in order to increase 
the effective area of the transient search at the expense 
of increased systematic errors in determining transient 
flux densities. 

3.5. Removal of Spurious Sources 

As discussed in Section |1.3| dividing up the data be- 
fore imaging can result in a decrease in the number of 
spurious sources (particularly those due to image defects 
that do not obey Gaussian statistics). This helps in the 
detection of transients, but also helps to generate purer 
single-epoch catalogs. We therefore accept only sources 
that are independently detected in both frequencies in 
at least one epoch (with a position matching tolerance of 
50", corresponding to a false match probability of < 2%). 
In many cases such sources are detected in additional 
epochs, although not necessarily always at both frequen- 
cies in every epoch, particularly if they are faint or if a 



single-epoch image is of poor quality. 

Our threshold of ~ 4.2a for detection in a single image 
corresponds to a threshold of ~ 5.9a in the dual-image 
catalog. We generated catalogs for the deep fields, con- 
sisting only of sources detected at both frequencies. The 
catalogs for the ELAIS-N1 (Table [2), Lockman (Tablet^, 
and Coma (Table E| deep fields contain 238, 189, and 186 
sources respectively. The new catalog for the NDWFS 
deep field (Tablej5j contains 195 sources. 

Examples of delects successfully removed from our 
data by the two-frequency requirement are spurious de- 
tections near bright sources, amplified noise close to the 
edges of images, and other failures of the source finder, 
including those due to underestimates of the background 
rms (e. g. , Fig.[9|. Although this method will not remove 
all spurious sources, examination of our data shows that 
it is successful enough to greatly reduce the number of 
false positive sources detected, which is extremely helpful 
in the search for transients. 

This method does not result in the removal of all spu- 
rious transient candidates, however. An example of a 
strong transient candidate seen in our data is shown in 



Fig. 10 — the candidate is present in images at both 3040 
and 3140 MHz. Further examination, however, showed 
that this apparent transient is an image of a bright source 
at a nearby position, and that the data were erroneously 
included in this image due to a telescope pointing error. 
Similar kinds of pernicious errors were also noted in VLA 
data by |Ofek et al. | ( [2010] ) , |Thyagarajan et al.| ( [20TT| , and 
|Frail et al.) ( |2U12| ), so even transient candidates that sur- 
vive all other cuts must be considered with appropriate 
skepticism until they are examined in detail. 

3.6. Completeness 

As in Paper II, we can assess the completeness and 
reliability of our data by plotting the source count his- 
tograms of the individual epochs and comparing to those 
from the master catalogs. The small size of the field com- 
pared to ATATS results in a larger fraction of sources de- 
tected in regions where the sensitivity falls off due to the 
primary beam. We therefore compute completeness con- 
sidering only sources from central regions of the mosaic 
with gain equal to one. 

For each epoch, we define completeness as 



C(S) = 1 



N U (S) 

N a (S) 



(1) 



where N U (S) is the number of sources in the master 
catalog (from regions with primary beam gain equal to 
one) with flux density brighter than 5, and no match 
in the single-epoch catalog within 150"; and N a (S) is 
the total number of sources in the master catalog with 
flux density brighter than S. This statistic does not rely 
on an accurate measure of the flux densities in individ- 
ual epochs; it simply determines the fraction of sources 
in the master catalog that are matched in a particular 
epoch, as a function of master epoch flux density. For a 
given epoch, completeness tends to increase as flux den- 
sity increases (Fig. 11). For most epochs, completeness 
crosses a threshold of 90% at some point, which we de- 
termine to be the completeness value for that epoch. For 
a handful of particularly bad images, completeness never 
crosses the 90% threshold (e.g. , Fig. [12]). 
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Figure 8. Left: Slopes of RANSAC fits to the ELAIS 3140 MHz single epoch data when matched to the ELAIS 3140 MHz master catalog. 
A slope, or post-imaging calibration (PIC) gain of 1 indicates epochs where flux densities in the individual epochs were unbiased relative 
to the deep image. Large excursions represent epochs with poor fits, typically due to calibration problems or other issues with the array. 
Periods where array performance was not optimal (for example, the excursions around MJD 2455560, i. e. , late December 2010) can be 
seen. Checking the ATA email archive, we find that it was snowing at the observatory at that time, which was probably the reason for the 
poor calibration in this particular case. Right: Histogram of the PIC gains from the left panel. Epochs with good data tend to have PIC 
gains between about 0.9 and 1.1. 



Table 2 

ELAIS Nl catalog 



PiGSS ID 


R. A. 


Decl. 


b-maj 


bmin 




S deep 


(mJy) a 




n de t h 


xl c 


vs/S d 




(hr) 


(deg) 


(arcsec) 


(arcsec) 


(deg) 


3040 


3140 


3040 


3140 


both 3040 3140 


3040 3140 


J160007p540202 


16.00209 


54.03392 


100.0 


100.0 


0.0 


9.31 ± 0.90 


8.38 ± 0.94 


1 










J160013p550244 


16.00393 


55.04559 


100.0 


100.0 


0.0 


5.89 ± 0.89 


7.25 ± 0.89 













J160027p543634 


16.00775 


54.60948 


100.0 


100.0 


0.0 


6.10 ± 0.67 


4.90 ± 0.57 





1 







J160034p542253 


16.00949 


54.38160 


181.1 


137.7 


69.6 


22.66 ± 0.98 


24.43 ± 0.78 


2 


1 







J160037p542734 


16.01066 


54.45914 


100.0 


100.0 


0.0 


3.28 ± 0.47 


4.06 ± 0.51 













J160043p551327 


16.01203 


55.22435 


330.1 


101.8 


79.3 


26.54 ± 1.49 


14.89 ± 1.76 





1 







J160121p545319 


16.02259 


54.88886 


100.0 


100.0 


0.0 


3.19 ± 0.45 


2.50 ± 0.63 













J160126p543931 


16.02405 


54.65880 


100.0 


100.0 


0.0 


4.18 ± 0.35 


4.00 ± 0.44 





1 







J160127p535640 


16.02427 


53.94520 


115.0 


106.4 


69.1 


24.11 ± 0.93 


25.15 ± 0.97 


60 


29 


10 0.20 3.15 


0.11 0.36 


J160139p545516 


16.02777 


54.92278 


116.9 


109.8 


177.8 


10.28 ± 0.54 


11.09 ± 0.66 


2 











Note. — Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table, which shows all 
sources detected at both frequencies in the deep image, is available. 

a Flux densities as measured from the deep images made from all good data at each of the two frequencies, 3040 and 3140 MHz. Quoted uncertainties 
are the sum in quadrature of the uncertainties on the fit and on the background measurement as reported by SFIND. 

b Number of single-epoch images with completeness better than 50 m Jy in which this source was detected, at each frequency (3040 and 3140 MHz). 
Number of well-behaved epochs (PIC gain within ±0.1 of the median value; Section |3.4| ) with completeness better than 50mJy in which the source was 
simultaneously detected in each of the single-frequency images. 

c Reduced chi-squared, for the hypothesis of no variability, computed using only well-behaved epochs where the source is detected at both frequencies in 
5 or more epochs. 

d Standard deviation divided by the mean of the single-epoch flux densities, computed using only well-behaved epochs where the source is detected at 
both frequencies in 5 or more epochs. 
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Table 3 

Lockman catalog 



PiGSS ID 


R. A. 


Decl. 


bmaj 


bmin 





Sdeep 


(mJy) a 




ndet h 




xl 


c 




/S d 




(hr) 


(deg) 


(arcsec) 


(arcsec) 


(deg) 


3040 


3140 


3040 


3140 


both 


3040 


3140 


3040 


3140 


J103814p583014 


10.63751 


58.50404 


104.0 


100.7 


94.0 


25.55 ± 1.20 


27.53 ± 1.00 


25 


13 


6 


0.80 


1.20 


0.15 


0.18 


J103834p580259 


10.64299 


58.04998 


100.0 


100.0 


0.0 


31.36 ± 1.24 


31.76 ± 1.10 


40 


35 


17 


0.92 


1.20 


0.17 


0.17 


J103856p575246 


10.64891 


57.87946 


141.9 


108.9 


145.9 


147.13 ± 2.78 


139.52 ± 2.33 


110 


107 


63 


8.33 


9.22 


0.21 


0.24 


J103900p590543 


10.65016 


59.09553 


100.0 


100.0 


0.0 


26.45 ± 2.13 


22.12 ± 2.04 


1 


1 


1 










J103911p581435 


10.65339 


58.24363 


100.0 


100.0 


0.0 


9.36 ± 0.69 


9.50 ± 0.77 


2 


3 


1 










J103913p591424 


10.65398 


59.24062 


100.0 


100.0 


0.0 


15.97 ± 1.72 


16.06 ± 2.28 



















J103922p582909 


10.65619 


58.48596 


100.0 


100.0 


0.0 


32.63 ± 0.81 


31.86 ± 0.72 


84 


69 


44 


1.85 


1.84 


0.19 


0.21 


J103942p580614 


10.66180 


58.10429 


104.2 


101.0 


69.1 


106.05 ± 1.45 


108.84 ± 1.47 


110 


110 


65 


2.85 


1.85 


0.09 


0.08 


J103953p583214 


10.66488 


58.53742 


100.0 


100.0 


0.0 


2.31 ± 0.52 


2.07 ±0.54 



















J103958p582350 


10.66636 


58.39744 


100.0 


100.0 


0.0 


5.18 ± 0.53 


4.03 ± 0.61 




















Note. — Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table, which shows all 
sources detected at both frequencies in the deep image, is available. 

a Flux densities as measured from the deep images made from all good data at each of the two frequencies, 3040 and 3140 MHz. Quoted uncertainties are 
the sum in quadrature of the uncertainties on the fit and on the background measurement as reported by SFIND. 

b Number of single-epoch images with completeness better than 50 mJy in whi ch this source was detected, at each frequency (3040 and 3140 MHz). Number 
of well-behaved epochs (PIC gain within ±0.1 of the median value; Section |3.4[ ) with completeness better than 50 m Jy in which the source was simultaneously 
detected in each of the single-frequency images. 

c Reduced chi-squared, for the hypothesis of no variability, computed using only well-behaved epochs where the source is detected at both frequencies in 5 
or more epochs. 

d Standard deviation divided by the mean of the single-epoch flux densities, computed using only well-behaved epochs where the source is detected at both 
frequencies in 5 or more epochs. 



Table 4 

Coma catalog 



PiGSS ID 


R. A. 


Decl. 


bmaj 


bmin 




Sdeep 


(mJy) a 




U d et h 




xl c 






(hr) 


(deg) 


(arcsec) 


(arcsec) 


(deg) 


3040 


3140 


3040 


3140 


both 


3040 


3140 


3040 3140 


J125137p280644 


12.86017 


28.11004 


119.9 


102.4 


169.6 


6.47 ± 1.84 


14.74 ± 2.07 

















J125150p275119 


12.86404 


27.85549 


100.0 


100.0 


0.0 


61.31 ± 1.85 


59.30 ± 2.05 


10 


7 


2 








J125218p273219 


12.87162 


27.53900 


100.0 


100.0 


0.0 


9.59 ± 1.36 


9.63 ± 1.63 

















J125224p273613 


12.87360 


27.60361 


100.0 


100.0 


0.0 


12.27 ± 1.89 


6.55 ± 1.96 

















J125252p282224 


12.88133 


28.37352 


100.0 


100.0 


0.0 


10.13 ± 1.39 


9.05 ± 1.57 


1 


1 











J125254p273141 


12.88200 


27.53374 


100.0 


100.0 


0.0 


5.67 ± 1.14 


5.78 ± 1.16 

















J125306p274607 


12.88522 


27.76888 


105.7 


101.2 


179.6 


97.19 ± 1.60 


95.19 ± 1.26 


37 


36 


23 


2.21 


2.30 


0.15 0.18 


J125314p273944 


12.88734 


27.66226 


105.4 


103.4 


171.0 


164.66 ± 2.43 


150.22 ± 1.42 


38 


37 


23 


12.88 


1.64 


0.18 0.09 


J125316p280328 


12.88842 


28.06191 


145.1 


106.9 


28.0 


5.63 ± 0.85 


7.98 ± 0.78 

















J125324p281827 


12.89001 


28.30763 


114.3 


101.6 


45.9 


9.60 ± 0.95 


7.42 ± 1.18 


















Note. — Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table, which shows all 
sources detected at both frequencies in the deep image, is available. 

a Flux densities as measured from the deep images made from all good data at each of the two frequencies, 3040 and 3140 MHz. Quoted uncertainties are 
the sum in quadrature of the uncertainties on the fit and on the background measurement as reported by SFIND. 

b Number of single-epoch images with completeness better than 50 m Jy in whi ch this source was detected, at each frequency (3040 and 3140 MHz). Number 
of well-behaved epochs (PIC gain within ±0.1 of the median value; Section |3.4| ) with completeness better than 50 m Jy in which the source was simultaneously 
detected in each of the single-frequency images. 

c Reduced chi-squared, for the hypothesis of no variability, computed using only well-behaved epochs where the source is detected at both frequencies in 5 
or more epochs. 

d Standard deviation divided by the mean of the single-epoch flux densities, computed using only well-behaved epochs where the source is detected at both 
frequencies in 5 or more epochs. 
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Table 5 

NDWFS catalog 



PiGSS ID 


R. A. 


Decl. 


h 

^TYl QjJ 


h 




q 

deep 


(mJy) a 




n d t° 




X 


2 c 


<J s 


/S d 




(hr) 


(deg) 


(arcsec) 


(arcsec) 


(deg) 


3040 


3140 


3040 


3140 


both 


3040 


V 3140 


3040 


3140 


J142444p341819 


14.41246 


34.30544 


105.5 


100.6 


123.3 


39.70 ± 1.71 


35.04 ± 2.41 


15 


4 


3 










J142448p341000 


14.41334 


34.16667 


100.0 


100.0 


0.0 


7.55 ± 1.42 


5.33 ± 1.50 



















J142515p345256 


14.42094 


34.88238 


112.0 


102.8 


76.1 


37.53 ± 2.43 


26.92 ± 2.75 


2 
















J142516p341558 


14.42156 


34.26741 


136.0 


100.5 


26.5 


9.33 ± 0.94 


11.47 ± 1.40 


1 
















J142522p340938 


14.42299 


34.16064 


100.0 


100.0 


0.0 


7.31 ± 0.91 


5.66 ± 1.26 





1 













J142540p345835 


14.42789 


34.97645 


124.0 


111.8 


66.3 


104.21 ± 3.21 


84.29 ± 2.90 


48 


11 


9 


1.29 


1.00 


0.17 


0.18 


J142543p335535 


14.42868 


33.92650 


100.0 


100.0 


0.0 


124.18 ± 2.80 


121.76 ± 2.90 


67 


67 


59 


2.40 


2.25 


0.10 


0.13 


J142607p340429 


14.43540 


34.07408 


100.0 


100.0 


0.0 


25.98 ± 0.85 


26.72 ± 0.83 


37 


25 


20 


1.56 


1.49 


0.19 


0.18 


J142609p333941 


14.43592 


33.66152 


100.0 


100.0 


0.0 


11.52 ± 1.50 


10.28 ± 1.74 


1 


1 













J142620p335126 


14.43910 


33.85727 


100.0 


100.0 


0.0 


5.57 ± 0.99 


4.57 ± 1.17 




















Note. — Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table, which shows all 
sources detected at both frequencies in the deep image, is available. 

a Flux densities as measured from the deep images made from all good data at each of the two frequencies, 3040 and 3140 MHz. Quoted uncertainties are 
the sum in quadrature of the uncertainties on the fit and on the background measurement as reported by SFIND. 

b Number of single-epoch images with completeness better than 50 mJy in whi ch this source was detected, at each frequency (3040 and 3140 MHz). Number 
of well-behaved epochs (PIC gain within ±0.1 of the median value; Section |3.4[ ) with completeness better than 50 m Jy in which the source was simultaneously 
detected in each of the single-frequency images. 

c Reduced chi-squared, for the hypothesis of no variability, computed using only well-behaved epochs where the source is detected at both frequencies in 5 
or more epochs. 

d Standard deviation divided by the mean of the single-epoch flux densities, computed using only well-behaved epochs where the source is detected at both 
frequencies in 5 or more epochs. 




Figure 9. Left: Part of the single-epoch 3040 MHz image for PiGSS Lockman 2010 April 7, with source positions from SFIND marked 
as circles of radius scaled with measured flux density. The image is 90% complete to llmJy. Right: Same region of the 3140 MHz image 
from the same epoch, also 90% complete to llmJy. A cluster of faint spurious sources is visible where the SFIND algorithm failed, and 
misidentified noise peaks as real sources. By keeping only sources that are present at both 3040 and 3140 MHz, we can exclude spurious 
sources in cases like this, and make our catalogs much cleaner. Examples like this one are rare in our data in the sense that they make 
up only a small fraction of detected sources, but if not identified and excluded, they dominate lists of transient candidates. Requiring 
two- frequency detections also excludes some faint real sources, such as the one at 10 h 50 m , 59° 11' in the left panel (also visible but not 
detected above the chosen threshold in the right panel), but in essence this is equivalent to increasing the overall threshold for detection. 
Greyscales are to lOmJy beam -1 . 



4. TRANSIENTS 

4.1. Comparing individual epochs to the deep image 

We searched for sources detected in one or more epochs 
(at both frequencies in at least one of these epochs) with- 
out a counterpart at either frequency in the deep image. 
As previously noted, the requirement of a detection in 
both frequencies dramatically lowers the false detection 
rate while preserving true transients. Nevertheless, any 
candidates must be examined carefully. 

4.1.1. Lockman 



In the Lockman field, a single transient candidate was 
detected, with flux density ~ 200 mJy, in nine (non- 
consecutive) epochs, but not in the catalog made from 
the deep image. Examination of postage stamps shows 
that the source is in fact clearly visible in all single epoch 
images, and the deep image, but was close to the edge 
of the mosaic, in a region where the mosaic gain cor- 
rection is a factor ~ 5. Small variations in the centroid 
position (and hence the gain) resulted in the source not 
meeting our cutoff of a factor 5 gain correction in many 
epochs (despite being detected by the source finder in 89 
of them), but just meeting it in the nine where it was 
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Figure 10. Subregion of the single-epoch 3040 MHz image for PiGSS ELAIS-N1 2010 October 5 (left). Comparison with the same 
frequency on epoch 2010 October 7 (right) shows a bright transient candidate (which is not present in any other epoch of ELAIS-N1, but 
is also present in the 3140 MHz image on 2010 October 5). More detailed examination shows that this is due to the inclusion of a small 
amount of data where the telescope was not pointed correctly. Greyscales are to 50mJy beam -1 . 



09-04-10 C90=11.335 




Figure 11. Left: Completeness, C(S) = 1 — (N u (S)/N a (S)), i. e. , the fraction of sources in the master catalog brighter than flux density 
S (in a region with gain equal to one) with a match in the single epoch catalog for ELAIS-N1 3140 MHz epoch 2010 September 4. The 
horizontal line shows 90% completeness, and the error bars were calculated following Camer on] ( |2011| ). The image quality in this field is 
typical (see Fig. [2J, and 90% completeness is reached at 12mJy. Right: Image of the field. The greyscale runs from to lOmJy beam -1 . 



detected. We conclude that there are no transients of 
one or more epochs in Lockman without a counterpart 
in the deep image. 

We also searched for sources with counterparts in the 
deep image that brightened in any epoch to 10 times or 
more the measured flux density in the deep image. As 
we argue in Section |1.2[ such sources ought to be con- 
sidered along with transients when calculating statistics. 
No such sources were found in the Lockman field. 

4.1.2. ELAIS-N1 

In the ELAIS-N1 field, 16 transient candidates were de- 
tected. Six were in a single epoch that showed imaging 
defects that were coincident in the 3040 and 3140 MHz 



data, and so are not astronomical. Five were in another 
epoch where the source finder failed near a bright source 
at both frequencies, and spurious sources coincided. Two 
of the candidates were in regions where the primary beam 
gain correction was a factor ~ 5. The remaining three 
candidates occurred due to inclusion of data from mis- 



pointed observations (e.g., Fig. 10). We conclude that 
there are no transients of one or more epochs in ELAIS- 
Nl without a counterpart in the deep image. 

Two sources in ELAIS-N1 appeared to brighten to 10 
times or more their flux density in the deep image. In 
one case a 30mJy source in the 3040 MHz deep image 
was matched to a nearby faint (4mJy) source in the 
3140 MHz deep image. The fainter source was not de- 
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09-27-10 C90=99999 




Figure 12. Left: Completeness, C(S) = 1 — (N U (S) / N a (S)) , i. e. , the fraction of sources in the master catalog brighter than flux density 
S (in a region with gain equal to one) with a match in the single epoch catalog for E LAIS-N1 3140 M Hz epoch 2010 September 27. The 
horizontal line shows 90% completeness, and the error bars were calculated following | Cameron] ( |201l) . The image quality in this field is 
unusually poor (see Fig. [2|, and 90% completeness is never reached. As a result this field was rejected from our analysis. Right: Image of 
the field, where the defects causing poor completeness are obvious. This epoch would be easily rejected upon manual inspection, but in 
large surveys where it is not possible to examine each image, completeness proves to be a useful method for rejecting bad images such as 
this one. It would be possible to implement an iterative approach to rejecting such data when constructing the deep image, but in practice 
the 3.5 mjy beam -1 rms cut proves adequate for this purpose. The greyscale runs from to lOmJy beam -1 . 



tected in the individual epoch 3140 MHz images, but the 
brighter source was (sometimes at > 40 mjy, since it was 
close to the edge of the mosaic where gain variations also 
play a part). This lead to the source being flagged as a 
candidate highly variable object, although in reality it is 
not. 

The second candidate was in fact two sources which 
were resolved in some epochs but not in others, leading 
to apparent large variability. We conclude that there are 
no sources in ELAIS-N1 which truly vary by a factor of 
10 or more. 

4.1.3. Coma 

In the Coma field, the transient search identified 19 
candidates without a match in the deep image. Of these, 
16 were present in the same single epoch. Many of the 
single-epoch images exhibited defects to varying degrees 
near the bright radio source Coma A, but this epoch was 
unusually bad, and 16 of the defects coincided at both 
frequencies and so were identified as transient candidates. 
Another two of the 19 candidates also appear near to 
Coma A in other epochs. On visual inspection, they ap- 
pear to be sidelobes (they appear in clusters of other ap- 
parent sidelobes, although these candidates are the only 
two that coincide in position at both frequencies). The 
remaining candidate is close to the edge of the mosaic 
where the primary beam gain correction is a factor ~ 5. 
Once again, none of these candidates survives further ex- 
amination. 

Three sources in Coma appeared to brighten by more 
than a factor 10 compared to their deep field flux densi- 
ties. Two of these candidates were close pairs that were 
resolved at some epochs but not in others; summing 
up the fluxes of the components, very little variability 
is in fact seen. The third candidate (where a lOmJy 
source appeared to brighten to > 100 mjy) occurred in 



the same bad epoch (with 16 transient candidates) men- 
tioned above, and is clearly an imaging defect related to 
Coma A. Once again, there are no highly variable sources 
in our Coma data. 

4.1.4. NDWFS 

In the reprocessed NDWFS field (with corrected as- 
trometry and a slightly different reduction pipeline in- 
cluding the choice of CLEAN iterations), surveying a 
smaller number of epochs than in Paper II (due to the 
completeness threshold we apply here) and the slightly 
different method of transient identification we use here 
(including a larger match radius), we find a single 11 mjy 
transient candidate. Further investigation shows that 
this source was detected at 3 mjy in the 3040 MHz deep 
image, and undetected, although just apparent on visual 
inspection, in the 3140 MHz deep image, and so it is not 
truly a transient, but may have varied by a factor ~ 4. 

One source in NDWFS appeared to brighten to 10 
times or more its flux density in the deep image. In this 
case a bright (300 mjy) source in the 3040 MHz deep im- 
age was matched to a nearby faint (llmJy) source in 
the 3140 MHz deep image. The fainter source was not 
detected in the individual epoch 3140 MHz images, but 
the 300 mjy source was, making it appear to be highly 
variable. Our chosen radius for source association is 
comparable to the synthesized beam size, so such mis- 
matches are not common in our data. However, cases 
like this could be better discriminated against by fit- 
ting Gaussians in every single-epoch image with posi- 
tions constrained to match the sources detected in the 
deep image. We intend to implement such a feature in 
our pipeline for future analysis. 

Once again, we find no highly variable sources in this 
field when comparing single epochs to the deep image. 
This result is in agreement with the lack of transients 
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seen when using a > 6<r threshold on the NDWFS anal- 
ysis presented in Paper II. 

4.2. Comparing monthly images to the deep image 

We searched for sources detected in one or more of the 
monthly images (at both frequencies in at least one of 
these epochs) without a counterpart at both frequencies 
in the deep image. 

In ELAIS-N1, seven candidates were found. Six were 
positionally coincident with the single epoch transients 
discussed above, which as previously noted, are not true 
astronomical transients. The seventh was a marginal de- 
tection at both frequencies in 2010 December, and at one 
frequency in 2010 October. Inspection of postage stamps 
shows that it was, in fact, also faintly visible in the other 
monthly images and the deep image. There are therefore 
no true transients in the ELAIS-N1 monthly data. 

In Lockman, two candidates were found. 
PiGSS J103818+584118 is visible only in the 2009 
December image, with 63040 = H-8 ± 2.1 mJy and 
63140 = 9.5 ± 2.1 mJy. The probability of such a source 
being spurious, assuming Gaussian statistics, is 2 x 10 -6 , 
although as previously noted, image defects and other 
issues mean that spurious sources do in fact arise more 
often than pure Gaussian statistics would imply, and in 
cases where the uncertainties are underestimated, the 
significance of a detection may be overstated. However, 
this candidate is detected with a joint significance of 
> 7a, which suggests that it is probably astronomical. 
The source is not extended relative to the restoring 
beam. The mosaic primary beam gain correction at 
the source position is a factor 2.8. The source is not 
obviously detected on visual inspection of postage 
stamps of individual epoch images from 2009 December, 
which is to be expected since its measured flux density 
(uncorrected for the primary beam) is ~ 4mJy, below 
the typical single-epoch completeness limit. Postage 



stamp images of the monthly data are shown in Fig. 13 



The closest NVSS source is over 6' away. The source 
apparently brightened by a factor > 6 relative to the 
2cr upper limit from the deep image (1.6 mJy beam -1 ), 
so is not a confirmed transient as we define them here, 
in the sense of a source that can be shown to have 
brightened by a factor of at least 10. The upper limit 
from FIRST is 1.02 mJy beam -1 , however, suggestive of 
a brightening of a factor > 10 if the source is assumed to 
have a flat spectrum. The candidate is 40" away from a 
Giant Meterwave Radio Telescope (GMRT) source seen 
at 610 MHz with a flux densi ty of 1.189 ± 0.137 mJy, 
GMRTLH J103820.1+584153 ( (Garn et al.][2008l ), but 
the areal density of GMRT sources is quite high, and we 
calculate a 30% probability that this association is due 
to chance. The closest Sloan Digital Sky Survey (SDSS) 
galaxy is also ~ 40" away from our transient candidate. 

The second candidate in the Lockman monthly data, 
PiGSS J104741+572316, appeared in 2010 February but 
is quite faint (63040 = 4.6 ± 1.1 mJy; 63140 = 4.8 ± 
1 . 1 m Jy ) . Visual inspection of the postage stamps 
(Fig. 14 ) suggests that it is probably not real. It appears 
similar to other background fluctuations in the image, 
particularly at 3140 MHz where image quality is notice- 
ably worse than for other months. 

In Coma, 13 candidates were found. One was position- 
ally coincident with a single-epoch transient discussed 



above, and so is spurious. Another 9 were found in a 
single image (2011 January) that showed a rather large 
number of presumably spurious sources of similar bright- 
ness (typically a few m Jy ) . Several had flux ratios in the 
two bands that would imply unrealistic spectral indices. 
All, except for the candidate that was coincident with the 
single-epoch candidate (which we also reject) do not ap- 
pear to be real candidates after further inspection. Once 
again the presence of the bright source Coma A in the 
field makes calibration difficult and leads to large num- 
bers of faint sidelobes or other imaging defects in the 
images, some of which are picked up by our transient 
pipeline. 

In NDWFS, seven candidates were found. All bar one 
were from regions of the image with primary beam gain 
> 2. On visual inspection, all appeared similar to other 
background fluctuations in the images. 

No sources in the Lockman, ELAIS-N1 or Coma fields 
were observed to vary by more than a factor of 10 in 
the monthly images compared to the deep images. One 
source in the NDWFS field was detected at 282 mJy in 
the deep 3040 MHz image, and matched to a nearby 
llmJy source in the 3140 MHz image. The source was 
detected at ~ 280 mJy at both frequencies in all NDWFS 
monthly images, but the match to the 11 mJy source re- 
sulted in it being flagged as a candidate highly variable 
source, although once again it is not truly highly variable. 

In summary, we find one transient candidate (although 
not a confirmed transient that is shown to vary by more 
than a factor 10) when comparing the monthly images 
to the deep images. 

4.3. Comparing the deep catalogs to NVSS 
We also compared the catalogs generated from the deep 



images (e.g., Tables 



J2||3||4| and§ to the NVSS cata- 
log. We used a matcn radius of 4o> , chosen to produce 
an expectation of less than one false match between the 
catalogs. 

With the corrected astrometry, and our more strin- 
gent detection requirements for the final catalog, 99% of 
PiGSS NDWFS sources have a match in NVSS. 

The flux limit of the NVSS is - 2.5 mJy, so flat 
spectrum PiGSS sources fainter than 2.5 mJy, or PiGSS 
sources fainter than ~ 1.4 mJy with typical spectral in- 
dices, a « —0.7 (Sv oc v a \ will not have NVSS counter- 
parts. We find three PiGSS sources > 2.5 mJy without 
NVSS counterparts in the Lockman field, five in ELAIS- 
Nl, five in Coma, and two in NDWFS. Only a total of 
five of these sources are brighter than 5mJy in PiGSS; 
the brightest candidate is llmJy. We examined postage 
stamps for all of the candidates, and all except two have 
faintly visible counterparts in the NVSS images, which 
were too faint to be detected in the NVSS catalogs. Of 
those two, one is clearly a sidelobe of Coma A. 

The remaining source, PiGSS J104711+582817, is seen 
with a flux density of 10.8 ± 0.5 mJy in both the 3040 
and 3140 MHz deep images of the Lockman field. It is 
also detected in over 20 single epoch images with consis- 
tent flux (Fig. 15), and visible on inspection in almost 
all epochs. The daily flux density measurements show 
no clear trend with time. There is no counterpart in the 
NVSS catalog (which is complete to 2.5 mJy), and noth- 
ing on visual inspection in the NVSS images. We mea- 
sured a 2cr upper limit of 2.5 mJy beam -1 at the ATA 
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Figure 13. Postage stamp images of Lockman monthly transient candidate PiGSS J103818+584118. The deep image is shown in the 
first panel, and then images from 2009 December, 2010 January, February, March, April, and May. Data are at 3040 MHz (top row) 
and 3140 MHz (bottom row). The greyscale runs from to lOmJy beam -1 . The transient candidate is visible in the center of the 2009 
December postage stamp at both frequencies. The upper limit from the deep i mage (< 1.6 mJy beam -1 ) is not faint enough to determine 
whether the source brightened by a factor > 10, to meet our definition (Section|1.2| of a transient. 
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position from the NVSS images, consistent with the lack 
of detection in the NVSS catalog. A search of the FIRST 
catalog shows a counterpart at 3.4 ± 0.15 mJy; although 
the source is present in FIRST and absen t from NVSS, 



Figure 14. Postage stamp images of Lockman monthly transient candidate PiGSS J104741+572316. The deep image is shown in the 
first panel, and then images from 2009 December, 2010 January, February, March, April, and May. Data are at 3040 MHz (top row) and 
3140 MHz (bottom row). The transient candidate is visible in the center of the 2010 February postage stamp at both frequencies. It appears 
similar to other background fluctuations in the image, and so the significance of the detection is likely overestimated, and it is probably 
not a real source. 

vary by more than a factor 10 between NVSS and PiGSS. 
If we assume a flat spectrum, however, and simply com- 
pare the flux densities from NVSS and PiGSS, we find 
no sources that varied by a factor 10 or more. 

Reversing the sense of the match, and considering only 
NVSS sources in the region where the PiGSS primary 
beam gain is unity, we find a total of 18 NVSS sources 
brighter than lOmJy without a cataloged PiGSS coun- 
terpart in the four deep fields. Of the 18, 10 are in the 
Coma field. On inspection, all in fact do have counter- 
parts in the PiGSS images, and the lack of matches is, in 
all cases, due to different deconvolution of source struc- 
tures at the two frequencies. Many of the NVSS sources 
fainter than lOmJy also appear to have counterparts in 
the PiGSS images, although some are truly undetected, 
due to the fainter flux limit of NVSS compared to PiGSS 
for sources with normal spectral indices. 

4.4. Rate limits 
Following |Bower et al. (2007), the surface density (or 



i t was not included in the analysis of |Levinson et al. 
(2002]) because it is fainter than the threshold they use 
to select t ransient cand i dates . The 325 MHz VLA obser- 
vations of Owen et al. ( 2009[ ) detected a counterpart at 
2.88 ± 0.16 mJy, a nd the 1)10 MHz GMRT observations of 
|Garn et aL (2008) detected a 2.28±0.14mJy counterpart. 
The radio sources are associ ated with a z = 0.55 QSO, 
SDSSJ104711. 15+582820.7 ( |Richards et la7|[2QQ9] ). The 
FIRST and NVSS flux densities suggest that the source 
exhibits variability. This would also explain the bright 
detection in PiGSS, although this may also be partially 
due to an inverted spectrum between 1.4 and 3.1GHz; 
since we do not have simultaneous radio detections at 
multiple widely-spaced frequencies we are unable to de- 
termine an instantaneous spectral index for the source. 

The degeneracy between variability and spectral index 
also makes it difficult for us to clearly detect sources that 



"two-epoch rate" ) tor a survey with N e epochs that cover 
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Figure 15. Left: Lightcurve for PiGSS J104711+582817, which appears to be a variable z = 0.55 AGN. 3040 MHz flux densities are shown 
in green, and 3140 MHz in blue. The deep field flux densities are shown as horizontal lines, with grey boxes representing the measured 
uncertainties on the deep field flux densities. Center: NVSS postage stamp. The green circle shows the position of the ATA source with 
no match in the NVSS catalog. The greyscale runs from to lOmJy beam -1 . Right: ATA 3040 MHz deep field postage stamp. The green 
circle shows the position of the ATA source with no match in the NVSS catalog. The greyscale runs from to 10 mJy beam -1 . 



a total area A to sensitivity S is 



S(>5) 



(N e - 1)A(>S) 



(2) 



where N t is the number of transients detected. Where 
no transie nt is detected , the 95% confidence upper limit 
is N t = 3 ( |Gehrels|1986| >. Here, A(>S) refers to the solid 
angle over which a source of flux density S or greater can 
be detected. Since faint sources are easier to detect close 
to the mosaic center than in regions where the primary 
beam gain falls off, we must compute A as a function of 
flux density. 

We use the 90% completeness values determined in the 
central regions of the mosaic where gain equals one to de- 
termine the detection limit for transients in a particular 
epoch. We then assume that the detection limit gets 
higher in inverse proportion to the gain falloff towards 
the mosaic edges, as the area above a certain primary 
beam gain increases. To determine transient limits for 
the survey as a whole, we sum the effective survey areas 
across epochs, as a function of sensitivity. 

For PiGSS sources missing in NVSS, corresponding to 
transient candidates on timescales of months to years, we 
set an upper limit to the transient rate (given that we 
see no transients brighter than 2.5 mJy) for long duration 
transients of lOmJy, T^i ong < 0.08 deg -2 . Rate limits for 



such sources are shown as the magenta curve in Fig. 16 



For PiGSS transients seen in the monthly data when 
compared to the deep image, of which we see one possible 
candidate, although no confirmed transients, we obtain 
^monthly ;$ 0.02 deg -2 at lOmJy. Rate limits for these 



sources are shown as the cyan curve in Fig. 16 



For PiGSS transients seen in the daily data when com- 
pared to the deep image, of which we see none, we obtain 
^daiiy ^ 0.009 deg -2 at lOmJy. Rate limits for these 
sources are shown as the blue curve in Fig. 16 



4.4.1. Rates per galaxy in Coma 

The Coma field contains a large over density of galax- 
ies, and is therefore more likely to contain supernovae 
than a random extragalactic field. At the distance 
of Coma (100 Mpc, corresponding to the cluster cen- 
ter at ~ 7000 km s~ ), our median 90% completeness 
in the single-epoch images, 15.2 mJy, corresponds to 
~ 4 x 10 28 erg s -1 Hz -1 . This is comparable to the peak 



radio luminosity of t he unusually luminous SN1998bw 



(Kulkarni et al. 1998). 



To estimate the number of galaxies in the Coma cluster 
surveyed by our observations, we searched the NASA Ex- 
tragalactic Database for galaxies with redshifts between 
5000 and 9000 km s _1 , within 60' of the center of our 
Coma map. This corresponds to approximately the area 
with mosaic gain unity in the deep image. We found 
665 galaxies within this area, which is consiste nt with 
the Coma clu ster galaxy density determined by |The &; 
Whitej ( [1986b - 

In the 73 days from the beginning to end of our Coma 
observations, no transients were seen. Taking a 2a upper 
limit of 15 transients per year, we can therefore constrain 
the rate of SN 1998bw like events to < 0.02 yr -1 per 
average galaxy in the Coma cluster. 

The average supernova rate in local galaxy clusters 
(including Coma) is a pproximately 0.01 galaxy -1 yr -1 
(Mannucci et al. 2008). The fraction of SNe detectable 
in the radio is still poorly constrained, and depends on 
survey sens itivity, but is typically greater than 10% (|Lien 



et al. 2011). This predicts that the radio transient rate 
from SNe in Coma ought to be > 0.001 galaxy -1 yr -1 . 
With order of magnitude or more improvements (com- 
pared to PiGSS) in sensitivity, area, or length of the 
observing campaig n, seeing ordinary radio supernovae 
(Weil er et al.||1998| in such a blind search might be pos- 
sible. 



5. VARIABILITY 

The PIC gain calibration (Section |3.4| ) helps us exclude 
apparent variability that is in fact due to bad calibration 
for particular epochs. For example, only one source in 
Lockman (Fig. [Tt]) was detected with a flux density that 
brightened to more than twice its flux density in the deep 
image. In this case the brightening was seen in only one 
of the two frequencies, which immediately arouses suspi- 
cion that it is not intrinsic to the source. In fact, the PIC 
gain fit failed to converge for this epoch and so the epoch 
was excluded from our variability analysis. The PIC gain 
cuts remove a number of similar cases where sources that 
appear highly variable in an individual epoch in fact have 
poor calibration for all sources in that epoch. 

Nonetheless, some erroneous flux densities still con- 
tribute to apparent large variability. Fig. [18] shows the 
source with the fourth highest variability in the Lockman 
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Figure 16. Cumulative two-epoch source surface d ensity for ra dio transients as a function of flux density, based on Fig. 9 of Bower et al.| 
(120071), Fig. 20 of Paper I, and updated in the light oflFrail et al.M2012fc. Bl ack lines show upper limits or measurements from the literature: 
from the s tudy of highly variable sources in the Lockman Hole ax 1.4 GHz jCarilh et al.|2 003 C03); from the VLA a rchival observatio ns of 
|Bower et al.j with a lyr timescale (labeled B07L) and a 2 month times cale (BU7S); Irom UKb follow-up observatio ns (Frail et al. 2003 F03); 
irom structure function analysis of NVSS and FIRST | Ofek et al. |2011| Oll), reinterpreted in the light of Frail et al. (2012}; from a 325 MHz 
VLA Lockman Hole archival transient search ( | Jaeger et al.|2 U12,J12), p lotted here as a 2 a upper limit; from a 5 (Jrlz (Jalactic Plane survey 
(Becker et al. 2010 B10); from an archival survey ol VLA calibrators (|Bell et al. 2011 B ell): from a VLA archival survey of the 3C 286 
field (Bower & Saul [2011| BS11); from comparison of NVSS and FIR ST (Gal- Yam et al. 2006 GY06); from the MOST archive (Ban nister| 
|et al.|20TT] h>all); and Irom drift-scanning surveys ( |Matsumura et ah 2009 MU9 ). (Jreen fines correspond to upper limits from ATATS: the 
upper straight line fro m the comparison of the deep tielcl to NVSS jurorf et al.||2010| , and the lower curve from analysis of the ATATS 
single-epoch images by |Croft et al.| |2011a| >. The magenta curve shows the upper limit on long-duration transients present in our PiGSS 
catalogs but missing in IN VSS. The cyan curve shows the upper limit on PiGSS transients from our monthly images when compared to 
the deep image. The blue curve shows the upper limit on PiGSS transients from daily images compared to the deep image. All of the 
PiGSS results are combined limits from all four PiGSS deep fields. The red lines show the predicted rates for the various source classes 
(tidal disrup tion events. Type II SNe, orphan gamma-ray burst afterglows, neutron star - neutron star mergers, and collapsar-like events) 
discussed by |Frail et aT1 ( |2012| ). 

field (for sources with detections in at least five individual 
epochs), as measured by the maximum absolute differ- 
ence in flux density between any two epochs. Although 
some of the outlying points are rejected due to PIC gains 
outside our range of acceptable values, remaining outly- 
ing flux density measurements are still responsible for 
this source being identified as highly variable. However, 
in many cases, the PIC gain selection does improve our 
ability to accurately measure variability. 

For the remainder of our variability analysis, we con- 
sider only data from "good" epochs with PIC gains 
wit hin ± 0.1 of the median PIC gain for the field (Sec- 
tion 3.4), that also meet our completeness cut (Table [lj. 



for all sources detected in these good epochs. In Fig. [19_ 
we show a histogram of the measured modulation frac- 
tions. We see that, as expected, the majority of measure- 
ments are consistent with variability less than a fraction 
of -0.1. 

We also compute the maximum daily fractional mod- 
ulation in each source, 



MAX 



(4) 



5.1. Fractional modulation 
We compute the daily fractional modulation, 

Si-S 
mi = — = — 



(3) 



for sources detected in at least 10 good epochs. In 
Fig. [20] we show examples of sources with large and small 
mmax- There appears to still be a tendency for outlier 
points to be responsible for some of the largest degrees 
of detected variability, as determined by variability seen 
at one frequency but not at the other. In Fig. [2l1we plo t 
rrimax as a function of S. As noted by Ofek etaL| (2011 ), 
fractional modulation and related variability measures 
have irregular statistical properties that depend on the 
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Figure 17. Lightcurve for PiGSS J105159+593240. 3040 MHz 
flux densities are shown in green, and 3140 MHz in blue. The 
deep field flux densities are shown as horizontal lines, with grey 
boxes representing the measured uncertainties. This source ap- 
pears to brighten in one epoch to more than twice its flux density 
in the deep Lockman image, although the variation is seen only at 
3040 MHz and not at 3140 MHz, and is therefore not intrinsic to 
the source. Some of the variability appears correlated between the 
two frequencies, whereas some does not. Some of this variability is 
nevertheless within the expectations from Gaussian unce rtain ties, 
although there is also a long tail of outlier points (Fig. \22\ . In 
the example shown here, the ~ 100 mJy outlier is from an epoch 
where the PI C ga in calibration at 3040 MHz did not converge to 
a fit (Section |3.4| , and so the bad epoch is in fact automatically 
excluded from our variability analysis. 
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Figure 18. Lightcurve for PiGS S J105141+591305, a z = 0.44 
SDSS QSO ( [Schneider eT!I1[2007| . 3040 MHz flux densities are 
shown in green, and 314UMHz m blue. Those within the accept- 
able range of PIC gain for Lockman (between 0.835 and 1.035) are 
shown as solid symbols, and those outside the acceptable range are 
shown as open symbols. The deep field flux densities are shown as 
horizontal lines, with grey boxes representing the measured uncer- 
tainties. This source appears to be highly variable when the max- 
imum absolute difference in flux density between any two epochs 
is calculated, and despite the PIC gain selection rejecting some 
outlying points (as well as a few which appear, at least from this 
lightcurve, to be acceptable) the largest outliers at 3040 MHz are 
not rejected in this particular case. In many cases, however, selec- 
tion based on PIC gain enables us to better measure real variability. 
The deep field flux densities are significantly higher than those in 
the individual epochs; this appears to be because the source ap- 
pears marginally resolved in the deep images. 
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Figure 19. Daily fractional modulation, rrti (computed for 
"good" epochs; see text) as a function of mean flux density, S, 
for sources detected in all four fields. 

number of epochs of observation, and so this statistic is 
not easily comparable from one survey to another. How- 
ever, it can still provide a measure of how often variabil- 
ity of some magnitude is seen in an individual survey. 

The upper envelope of the points in Fig. [21] shows that 
variable sources in our dataset do not typically vary by 
more than a factor ~ 2 from epoch to epoch, and in the 
few cases where more extreme variability appears to be 
detected, further examination of light curves and postage 
stamp images shows that it is not intrinsic to the source 
(for example, because variability is seen at only one of 
the two frequencies). This suggests that our threshold 
of a factor 10 (Section |1.2[ ) for "true" transients (typi- 
cally explosive or impulsive events), as opposed to the 
upper envelope of "normal" variability, is reasonable. Of 
course, extreme variables are still interesting in their own 
right, and may represent sudden changes in accretion rate 
onto AGNs, for example. 

5.2. Consistency of variability in the two frequency 

bands 

Where similar patterns of variability are seen at both 
frequencies simultaneously, these are more likely to be in- 
trinsic to the source than patterns where the variability 
is not correlated at the two frequencies. However, varia- 
tions in overall gain calibration of the telescope, as well 
as other systematic effects, can also lead to variations in 
both frequencies simultaneously, so this property alone 
cannot be used to distinguish intrinsic variability from 
that due to calibration uncertainties. 

To examine whether differences between measured flux 
densities at the two frequencies are consistent with the 
measurement uncertainties, we compute the difference 
between the two flux densities as a multiple of the total 
uncertainty. We define this difference, 



Si, 3040 — Si, 3140 



V / (A^,3040) 2 + (AS i|3 l40; 



2' 



(5) 
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Figure 20. Lightcurve for PiGSS J105703+584717 (top), a source 
with large maximum daily fractional modulation, m macc , and 
PiGSS J104719+582116 (bottom^ a source w ith small Tfimax as- 
sociated with a z = 1.22 QSO (Schneider et al. 2007). 3040MHz 
flux densities are shown in green, and 314UMJlz m blue. The deep 
field flux densities are shown as horizontal lines, with grey boxes 
representing the measured uncertainties. 

for each good epoch, z, where a source is detected in 
both frequencies. The uncertainties at each frequency are 
the reported uncertainties in the fit and the background 
for each measurement, added in quadrature. 

In Fig. [22j we plot a histogram of Di. If the differ- 
ences in measured flux densities obey Gaussian statistics, 
this should match the standard normal distribution. We 
see that the agreement is quite good, although there is 
a tail made up of a few points which have larger dis- 
crepancies than would be expected from Gaussian statis- 
tics, and the data histogram is skewed slightly to the 
right, implying that sources are measured to be slightly 
brighter at 3040 MHz than at 3140 MHz. In fact, if we di- 
vide the 3040 MHz fluxes by the corresponding 3140 MHz 
fluxes for the 8430 dual-frequency detections of sources in 
the epochs considered in our variability analysis, we find 
that the sources are 1.9% brighter at 3040 MHz than at 
3140 MHz. If we assume no systematic errors in over- 
all calibration, this implies a mean spectral index of 
a = —0.6, consistent with expectations. 

5.3. Standard deviation 
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Figure 21. Maximum daily fractional modulation, m max as a 
function of mean flux density, S, for sources with 10 or more de- 
tections (in good epochs) in all four fields (using data for both 
frequencies). 
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Figure 22. Difference between flux density measurements at the 
two frequencies, as a multiple of measured uncertainties on the flux 
densities (solid histogram). The dotted line shows the standard 
normal distribution, scaled to match the sample and bin sizes of 
the solid histogram. 



We also compute the standard deviation divided by 
the mean flux density for all sources with 10 or more de- 
tections in good epochs, which we plot in Fig. [23] Once 
again, the sources showing the most extreme variabil- 
ity (for example, the sources with as/S = 0.71 and 
0.53 in this plot) are often due to outliers (in each of 
these two cases, a single deviant 3040 MHz point which 
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is not matched by a correspondingly deviant 3140 MHz 
point). However, many other sources which would oth- 
erwise show up as extreme variables in such plots benefit 
greatly from the PIC rejection in having their variabil- 
ity accurately measured. For example, the most extreme 
source using this criterion without PIC rejection would 
be PiGSS J160833p543457 with_<7s/S = 0.96; after we 
reject the bad epochs, its (Js/S drops dramatically to 
0.16. The median value of as/S is ~ 6% lower than it 
would be for the ensemble of all sources if we did not 
apply the PIC rejection. 

In Fig. 24 , we show as/S for the monthly data. Again, 
many of the outliers show apparent variability in only 
one of the two frequencies. Some apparently extremely 
variable sources are resolved into multiple components 
(or associated with confusing sources) at some epochs. 
Some appear variable due to the smaller number of de- 
tections we require here. 

5.4. Chi- squared 

We also compute a chi-squared statistic, given the hy- 
pothesis that the source flux density does not change 
from epoch to epoch, 



X 



Sj-S 

AS,; 



(6) 



For a source detected in n epoc h good epochs, we can 
compute a reduced chi-squared, xt = X 2 /( n epoch — 1), 
enabling us to compare the significance of the detection 
of variability for sources detected in different numbers of 
epochs. We compute xl independently for the 3040 MHz 
and 3140 MHz light curves for each source, considering 
only epochs where the source is detected at both fre- 
quencies (from the subsample of epochs that passed the 
criteria for our variability analysis). Fig. [25] compares 
at the two frequencies. For sources detected in a 
large number of epochs (which naturally tend to be the 
brighter sources), the two measurements of xl tend to 
agree. These bright sources do not necessarily have larger 
fractional variability; the x 2 is larger because we can de- 
tect variability more easily for these well-characterized 
sources. 

Some fainter sources, but also a handful of brighter 
sources, have discrepant measures of xi, sometimes ap- 
pearing variable at one frequency but not at the other. 
Again, this is typically due to erroneous outlying flux 
density measurements in small numbers of epochs. If 
we had considered only a single frequency, however, we 
might have concluded that large values of xl generally 
indicate significant variability. We appear, t houg h, to be 
in a similar regime as for transients (Section 1.3), where 



the non- Gaussian tail of outliers is in danger of domi- 
nating truly interesting sources showing large intrinsic 
variability. 

In Fig. 26, we plot xl against mean as/S for the 



same sources. The variability seen in faint sources is 
less well-characterized than for their brighter counter- 
parts (because they are less well detected) and so they 
tend to have lower values of xl- Faint sources with 
larger measured variability (as/S) have higher xl be- 
cause such variability becomes increasingly large relative 
to the measurement uncertainties on the flux densities. 



Bright sources tend to have larger xl on average, but a 
wide range in both parameters, suggesting that some are 
well-characterized (as measured by xl) as relatively non- 
variable, and others well-characterized as more variable. 

5.5. Structure Function 

We calculate s t ructu re functions for our sources follow- 
ing jLoveireFar] ([2008]) , 



(7) 



3,k 



where S u j and S Ui k are pairs of flux density measure- 
ments normalized by the mean flux density of the source 
over all epochs (treating the two frequency bands, v, in- 
dependently), and N T is the number of pairs of measure- 
ments with time lag r in bins equally spaced in log(r). 
Most of the structure functions for sources with dual- 
frequency detections in at least 10 good epochs appear 
quite flat, with normalized amplitudes typi cally between 
0.01 a nd 0.1, comparable to the results of |Lovell et al. 
( |2008[) on shorter timescales. 

In Fig. 27 we show structure functions, computed inde- 
pendently in our two frequency bands, for the 12 sources 
with xl > 5 in both bands, i. e. , those where variability is 
best detected (although not necessarily the sources with 
the largest measured variability; see Fig. 26). Again, 
most of the structure functions appear quite flat, al- 
though a few of the sources show signs of an upturn on 
long time lags, and a few also still show disagreement be- 
tween the two frequency bands, indicative of remaining 
bad data, despite us using only data where sources are 
detected at both frequencies in at least 10 good epochs 
(defined as PIC gain within the good range) here. The 
fiat structure functions likely indicate that intrinsic vari- 
ability of these sources is rather low, although some may 
be variable on timescales shorter than measured here (for 
example, as intraday variables). 

6. CONCLUSIONS 

We report on an ATA 3.1 GHz program that observed 
one of four 12 square degree deep fields an average of a 
little less than once per day, as part of the PiGSS project 
(covering all four fields over the course of ~ 2yrs). We 
detect no transients in comparisons of images made from 
379 individual epochs to deep images made from all us- 
able observations of a particular field, or from compar- 
ison of the deep images to images made from monthly 
averages of our data or to legacy surveys from the liter- 
ature. We place upper limits on the surface densities of 
transients on timescales of days to years, and flux den- 
sities brighter than a few mJy, that are still around an 
order of magnitude higher than predictions of rates from 
the l iterature, fo r exam ple for tidal disruption events 



from Frail et al.| ([2012]). Blind surveys covering more 
than an order ol magnitude greater effective area ought 
to begin to detect such events if theoretical predictions 
are correct. 

We discuss a number of caveats to be aware of in the 
search for transient and variable sources, and note in 
particular that the rejection of data from bad epochs 
(for example, by examining post-imaging calibration fits 
of the single-epoch flux densities to those from the deep 
comparison image, and by rejecting epochs with poor 
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Figure 23. Left: Standard deviation of the daily flux densities, as, divided by the mean flux density, S, for sources in_all four fields with 
10 or more detections in good epochs (using data from both frequencies), as a function of S. Right: Histogram of crs/S. 
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Figure 24. Left: Standard deviation of the monthly flux densities, as, divided by the mean flux density, S, for sources in all four fields 
with 3 or more detections in the monthly images, as a function of S (using data for both frequencies). Right: Histogram of crs/S. 



completeness) as well as rejection of image defects and 
other spurious sources which do not obey simple Gaus- 
sian statistics (in particular, by considering only sources 
that are detected in images made from neighboring fre- 
quency bands or from data split in some other way) are 
critically important in the regime where the potential 
false positive rate is large compared to the true rate of 
candidates of interest. We also note that some of the 
sources which appear most variable in our survey (and 
perhaps in others) appear variable because of discrepant 
flux density measurements in only one of our two frequen- 
cies, further reinforcing the dangers of relying on single 
flux density measurements without a simultaneous sec- 
ond confirming measurement. 



A forthcoming paper will discuss the results from a 
two-epoch survey, also undertaken as part of PiGSS, 
which surveyed thousands of square degree to similar sen- 
sitivity to the single-epoch images presented here. 
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Figure 25. Reduced chi-squared, x^, for the hypothesis that sources do not vary in individual epochs, computed for sources in (from 
top left to bottom right panel) the NDWFS, Lockman, ELAIS-N1 and Coma fields. We consider only epochs where sources are detected 
at both frequencies, and only sources with dual-frequency detections in at least 10 good epochs. Sources are color-coded according to the 
number of epochs in which they were detected. Circle sizes scale with mean flux density, S. The grey line guides the eye for sources with 
equal xt a ^ both frequencies. 
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Figure 27. Structure functions (see text) for the twelve sources where variability is best detected (xl > 5 in both bands). Data at 
3040 MHz (green) and 3140 MHz (blue) are considered independently. (Figure continued on next page . . . ). 
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